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In this study, the influence of intratumoral morphological heterogeneity of breast cancer on neoadjuvant 
chemotherapy (NAC) efficiency was investigated. In particular, we analysed the association of NAC 
response and pre- and post-NAC expression of the main multidrug resistance (MDR) genes - ABCB1, 
ABCC1, ABCC5, ABCG1, and ABCG2, with the presence of different morphological structures in breast 
tumors. In addition, the expression of MDR genes was investigated in different morphological structures 
and in their microenvironment by comparing probes obtained using laser microdissection. The results of 
this study showed that tumors with alveolar structures were more frequently N AC-nonresponsive than cases 
without this structural type (p = 0.0028, Bonferroni-corrected p = 0.014). The presence of trabecular 
structures in breast tumors was also associated with chemoresistance (p = 0.0272, Bonferroni-corrected p = 
0.136). High expression of MDR genes was not found in alveolar structures (including their 
microenvironment) and in tumors containing this structural type. In contrast, more active MDR genes and 
expression of the ABCB1 gene were found only in trabecular structures. Taken together, our data indicate 
that breast tumors with alveolar structures possess resistance to NAC, which is not related to high expression 
of MDR genes, whereas chemoresistance of tumors with trabecular structures can depend on the expression 
level of ABCB1. 



Breast cancer is a complex disease with high inter- and intratumoral heterogeneity 1,2 . Invasive carcinoma of 
no special type (IC NST) 3 , previously classified as invasive ductal carcinoma, not otherwise specified 
(NOS) 4 , account for a substantial proportion (up to 75%) of breast cancer cases and display extremely 
diverse morphological characteristics that make these tumors difficult to classify histologically 5 . IC NST tumors 
often contain minor components of special types of histology, including architectural features of the invasive 
process, which may vary widely both within a single tumor and from case to case 6 . Previously, we have described 
five different types of invasive component or morphological structures in IC NST tumors - tubular, trabecular, 
solid, alveolar structures, and discrete groups of tumor cells 7 - and found that such intratumoral morphological 
heterogeneity was related to cancer metastasis 7,8 . In addition, our recent data provided evidence that phenotypic 
drift can be cause of the development of intratumoral morphological heterogeneity in IC NST 9 . 

At present, the study of intratumoral heterogeneity appears to be key for the development of personalized 
approaches in cancer treatment 10,11 . In 2008, we reported that intratumoral morphological heterogeneity of IC 
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All patients had invasive carcinoma of no special type. 

Abbreviations: NAC, neoadjuvant chemotherapy; CAX, Cyclophosphamide-Adriamycin-Xeloda; 
FAC, 5-Fluorou rail-Ad riamycin-Cyclophosphamide; HER2 testing is performed in accordance with 
American Society of Clinical Oncology/College of American Pathologists Guideline 2007 
Recommendation 48 . 



NST was associated with the efficiency of neoadjuvant chemotherapy 
(NAC): the presence of alveolar and trabecular structures in breast 
tumors resulted in a poor response to NAC 12 . The mechanisms 
involved in drug resistance of heterogeneous tumors are not com- 
pletely established. It is reasonable to assume that tumors contain 
different clones of tumor cells with different degrees of responsive- 
ness to chemotherapy. Various molecular factors can be involved in 
tumor drug resistance, mainly ATP-binding cassette (ABC) trans- 



porters, which are encoded by a large family of multidrug resistance 
(MDR) genes and play a major role in mediating drug resistance 13 . In 
this study, we focused on the key MDR genes - ABCB1, ABCC1, 
ABCC5, ABCG1, and ABCG2. 

Thus, on a larger scale than in our previous research 12 , we aimed to 
study the association between the NAC response and the presence of 
different types of morphological structures in IC NST tumors. Then, 
we investigated pre- and post-NAC expression levels of the MDR 
genes in breast tumors and whether they depend on the presence of 
various morphological structures. Finally, using laser microdissec- 
tion we estimated expression levels of the MDR genes directly in 
different morphological structures and in their microenvironment. 

Results 

Analysis of correlation of intratumoral morphological hetero- 
geneity of IC NST and NAC response. Using hematoxylin & 
eosin staining and morphological analysis, we identified different 
morphological structures in breast tumors (n = 382; patient 
characteristics are listed in Table 1). Statistical analysis of the 
association between NAC response and the presence/absence of 
different morphological structures in breast tumors was conducted 
using Pearson's chi-square test with Bonferroni correction (Table 2). 
Patients with alveolar structures were more frequently NAC- 
nonresponsive than were cases without this structural type (61.9% 
vs. 46.4%; p = 0.0028, Bonferroni-corrected p = 0.014). In addition, 
breast tumors with trabecular structures more often demonstrated 
chemotherapy resistance as compared with tumors without these 
structures (58.8% vs. 45.3%; p = 0.0272). However, the difference 
did not reach statistically significance after Bonferroni correction (p 
= 0.136). 

Analysis of correlation of intratumoral morphological hetero- 
geneity of IC NST and expression of MDR genes. Using RT- 
PCR, we have compared the expression levels of MDR genes in the 
pre- and post-NAC tumor samples (n = 69) with the presence/ 
absence of different morphological structures. Statistical analysis of 
the correlation of gene expression with different morphological 
structures of tumor was performed using logistic regression with 
Bonferroni correction. Expression of ABCB1 (p = 0.007) and AB- 
CC5 (p = 0.027) genes was significantly lower in post-NAC tumors 
with tubular structures compared with tumors without tubular 
structures. Patients with solid structures displayed a decreased post- 
NAC expression of ABCBi (p = 0.011) and ABCG2 (p = 0.002) 
genes in comparison with patients whose tumors did not contain 
the solid structures. However, the differences in the expression of 
ABCB1 and ABCG2 were not significant after Bonferroni correction 
(p > 0.05; Table 3). In Tables 4 and 5, we showed the results of 
analysis of MDR gene expression in the distinct morphological 



Table 2 | The response to neoadjuvant chemotherapy dependin 
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Statistical analysis: p value, Pearson's chi-squared test; 

*, Bonferroni-corrected p value was calculated as the each p value multiplied by the number of tests (n = 5). 

Abbreviations: n, number of patients with presence (yes) or absence (no) of any morphological structures; N, number of patients with presenc< 
response/non-response to neoadjuvant chemotherapy; PR, partial response; SD, stable disease; PD, progressive disease; NA, not applied. 
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structures and in their microenvironment of two breast tumors 
obtained by laser microdissection. Tissue sections before and after 
laser microdissection are shown in Figure 1. In trabecular structures, 
we detected the expression of the most MDR genes. Most impor- 
tantly, expression of ABCB1 gene was specific only for trabecular 
structures. In contrast, alveolar structures displayed a low expres- 
sion of MDR genes. Interestingly, the microenvironment of different 
morphological structures also demonstrated expression of MDR 
genes, although in most cases gene activity was less than in morpho- 
logical structures. 

Analysis of correlation of intratumoral morphological hetero- 
geneity of IC NST and up-/downregulation of MDR genes 
during NAC. We performed the quantitative analysis of the MDR 
gene expression changes during NAC in breast tumors with different 
morphological structures. For this purpose, frequencies of the 
increase or decrease in MDR gene expression during NAC were 
compared between patients with the presence or absence of 
different morphological structures using Pearson's chi-square test 
with Bonferroni correction (Table 6). A decrease in ABCB1 
expression was more frequently detected in breast tumors with 
solid structures in comparison with tumors without this structural 
type (63% versus 35%, p = 0.022). Similarly, we frequently observed 
ABCG2 downregulation in cases with discrete groups of tumor cells 
than in patients whose tumors were free from these morphological 
structures (66% versus 32%, p = 0.008). However, the statistical 
significance of the differences was not confirmed after Bonferroni 
correction (p > 0.05). 

Discussion 

In this study, we demonstrated that intratumoral morphological 
heterogeneity of breast cancer, which was previously described in 
the most common histological type - invasive carcinoma of no spe- 
cial type 7 , influences the efficiency of neoadjuvant chemotherapy. On 
a larger scale than our previous study 12 , we confirmed, with greater 
statistical significance, the association between the presence of alve- 
olar structures in breast tumors and poor NAC response, which was 
significant after Bonferroni correction. In addition, trabecular struc- 
tures were also related to chemotherapy resistance as previously 
suggested 12 ; however, the Bonferroni correction did not showed 
the significance of this association. 

Chemotherapy efficiency is composed of many host factors and 
tumor alterations 14 . The most common reason for cancer drug res- 
istance is the expression of one or more ATP -binding cassette trans- 
porters that detect and eject anticancer drugs from tumor cells 13,14 . 
The data of our previous study suggest that changes in expression of 
MDR genes during the chemotherapy process or the development of 
adaptive MDR, but not the mRNA levels of these genes per se, are 
associated with NAC efficiency. In particular, reduction in MDR 
gene expression in post-NAC samples in comparison with pre- 
NAC tumors was linked with good response to chemotherapy, 
whereas patients displaying MDR gene upregulation exhibited res- 
istance to therapy 15 . In addition, recent data indicate that chemo- 
therapy-induced upregulation of MDR genes can result in decreased 
distant metastasis- 16 and disease-free 17 survival. 

The data obtained in this study suggested that chemoresistance of 
breast tumors with trabecular structures is associated with high levels 
of expression of MDR genes. Trabecular structures showed more 
active MDR genes. In addition, only trabecular structures were found 
to express the ABCB1 gene encoding P-glycoprotein, which is broad- 
spectrum drug efflux pump and plays a central role in MDR 14 . This is 
in consistence with our published study performed on the basis of 
FFPE tumor material 18 . Interestingly, an increased expression of 
MDR genes both before and after NAC, and their upregulation dur- 
ing chemotherapy were not found in breast tumors with trabecular 
structures. 
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Table 4 


Expression levels of MDR 
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21.594 


ABCG1 
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0 
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0.0004 
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0 


ABCG2 
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0.002 


0 
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0 


301.593 


0 


Abbreviations: ME, microenvironment. 



In contrast, poor response of breast tumors with alveolar struc- 
tures to NAC seems not to correlate to MDR gene activity. In par- 
ticular, tumors with alveolar structures did not show high expression 
of MDR genes both before and after NAC. Moreover, MDR gene 
upregulation was also absent in these tumors. In addition, laser 
microdissection-based expression analysis showed a low activity of 
MDR genes in alveolar structures. Previously, MDR genes and ABC 
transporters were found to be expressed in different cells of the tumor 
stroma 19 21 . However, we showed almost complete absence of MDR 
gene activity in the microenvironment of alveolar structures. 

Alveolar structures have spheroid shapes containing 10-30 tumor 
cells. During the 1970s, it was shown in vitro that multicellular 
spheroids containing thousands of cells were more resistant to radi- 
ation and adriamycin than cells cultured in a monolayer 22 ' 23 . 
Additionally, smaller spheroids consisting of 5-50 cells were also 
found to be more chemoresistant than single cells or monolayer 
cultures 24,25 . This phenomenon was later termed "multicellular res- 
istance (MCR)" 26 , which was linked to insensitivity to almost all 
anticancer drugs 27 . Subsequent studies reported that MCR mechan- 
isms are associated with the inability of the drug to penetrate the 
spheroid, differences in cell cycle distribution within spheroid, the 
presence of hypoxic cells, apoptosis inhibition, an increased DNA 
repair capacity, the so-called "contact effect", and others 25 ' 27,28 . In 
addition, formation of tumor spheroids is a hallmark in vitro of 
cancer stem cells, and these tumor- initiating cells demonstrate 
increased drug resistance and prolonged survival 25 . 

P-glycoprotein acts more efficiently in tumor spheroids than in 
monolayer cultures because of upregulation of the ABCB1 gene by 
hypoxia and acidification, which are attributes of tumor spheres 30 . 
For instance, Doublier et al. showed that tumor spheroids obtaining 
from MCF7 breast cancer cells were resistant to doxorubicin, and this 
resistance was associated with an increased P-glycoprotein express- 
ion via activation of hypoxia- inducible factor-1 (HIF-1) 31 . Based on 
the similarity between MCF7 spheroids and clusters of tumor cells of 
invasive micropapillary carcinoma (IMPC) of the breast, it has been 
hypothesized that chemoresistance of this subtype of breast cancer is 
provoked by high activity of P-glycoprotein 31 . IMPC clusters and 
alveolar structures of IC NST seem to have similar morphology. 
Moreover, micropapillary foci can be found in IC NST 32 . Results of 
our study demonstrate that high expression and upregulation of 
ABCB1 gene are not common for alveolar structures and breast 
tumors containing this morphological variant. The differences in 



chemoresistance mechanisms of IMPC clusters and IC NST alveolar 
structures probably reflect specific but not identified yet biological 
features of these groups of tumor cells. IMPC clusters detach from 
the stroma and display aberrant localization of glycoprotein MUC-1 
(also known as epithelial membrane antigen, EMA) at the stromal- 
basal surface, corresponding to an inversion of cell polarity 31,33 . In 
contrast, immunohistochemical analysis of two IC NST cases used in 
laser microdissection of the present study showed that in alveolar 
structures expression of MUC-1 (EMA) is observed in the whole 
cytoplasmic membrane and/or cytoplasm. In addition, in alveolar 
structures the change of cell polarization is not found, although the 
loss of interaction with stroma sometimes occurs (Fig. 2). It should be 
also pointed out that in IC NST, the frequency of alveolar structures 
comprises 60-83% (our unpublished data), whereas the proportion 
of micropapillary foci constitutes only 7.0% 32 . 

It is interesting to note that tumor spheroids (microemboli) were 
previously detected in peripheral blood of patients with lung 34-36 , 
prostate 36,37 , renal cell 38 , colorectal 39,40 , and breast 41 carcinoma. In 
comparison with circulating single tumor cells, such tumor clusters 
were shown to have anoikis suppression and the highest metastatic 
potential 35 . In addition, it has been suggested that the lack of apop- 
tosis and perhaps proliferation make tumor clusters more resistant to 
chemotherapy than solitary tumor cells 42 . 

Overall, our data confirm that intratumoral morphological het- 
erogeneity of invasive carcinoma of no special type is related to NAC 
efficiency. Breast tumors containing alveolar structures demonstrate 
a poor response to NAC and such observation is not explained by 
initial and adaptive MDR or upregulation of MDR genes during 
chemotherapy. In addition, the presence of trabecular structures in 
breast tumors is also associated with chemoresistance probably via 
ABCB1 expression (P-glycoprotein). Further studies are needed to 
investigate what factors/mechanisms are involved in chemoresis- 
tance of breast tumors with trabecular and alveolar structures. 

Methods 

Patients, tumors, and treatments. Patients (n = 382) with clinical stage IIA to IIIC 
(T^No^Mo) IC NST, between 25 and 71 years of age (mean age: 51.9 ± 0.50), and 
treated in the Cancer Research Institute (Tomsk, Russia) between 2006and2012 were 
included (Table 1). The procedures followed in this study were in accordance with the 
Helsinki Declaration (1964, amended in 1975 and 1983). This study was approved by 
the institutional review board, and all patients signed an informed consent for 
voluntary participation. All patients received two to four preoperative cycles of FAC 
(5-Fluorourail, Adriamycin, and Cyclophosphamide), CAX (Cyclophosphamide, 
Adriamycin, Xeloda) regimen, or Taxotere. Physical examination was performed 



Table 5 | Expression levels of MDR genes in different types of morphological structures and in their microenvironment of the second breast 
tumor 

Trabecular Discrete groups 

Genes Alveolar structures ME structures ME Solid structures ME of tumor cells ME 



ABCB1 0 0 2.210 0 0 0 0 0 

ABCO 00 00 000 0 

ABCC5 0.157 0 0 0.068 1260.523 0.543 0 0.181 

ABCG1 0.465 0.003 0.066 0.001 3.437 0.038 0.739 0.053 

ABCG2 0 0 21.360 0 0 0 29.912 0 
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Figure 1 | Laser microdissection of different types of morphological 
structures from sections of breast tumor, (a), alveolar structure; (b)> 
trabecular structures; (c), discrete groups of tumor cells (marked by 
arrows); (d), solid structures. 200X magnification. 

before NAC and was repeated after 2 cycles of NAC and before surgery to determine 
clinical response. Imaging of the primary breast lesion was performed with 
mammography and/or ultrasonography. Clinical and imaging responses to NAC 
were categorized into the following groups according to the International Union 
Against Cancer criteria 43 : complete response (CR), partial response (PR), stable 
disease (SD), and progressive disease (PD). Patients were grouped into clinical 
responders (PR) and non-responders (SD and PD). Fresh tumor tissues (n = 382) 
obtained after NAC were fixed in 10% neutral formalin (Karbolit, Russia) for 24 
hours, rinsed with an isopropanol solution (Biovitrum, Russia), and embedded in 
paraffin (Biovitrum, Russia). In addition, out of 382 fresh tumor samples, 69 
(randomly selected) samples and their biopsies (before NAC) were placed in RNA- 
later solution (Ambion, USA) and were stored at -80 "C until RNA isolation. The 
operative samples from two patients without NAC were collected after surgery in 
nitrogen and were stored at -80°C until laser microdissection. Cases with CR were 
excluded from this study because of a loss of post-NAC tumor samples for expression 
analysis. 

Morphological analysis. Morphological analysis included identification of different 
morphological structures in breast tumors and was performed by light microscope 
(Axio Scope, Carl Zeiss, Germany). The presence of tubular, trabecular, solid, alveolar 
structures, and discrete groups of tumor cells was evaluated in breast tumors (n — 
382) using 5 um-thick tumor sections stained by hematoxylin (Dako, Denmark) and 
eosin (Dako, Denmark). All tumor slides (five from each sample) were reviewed by 
three experienced pathologists. Tubular structures used in tumor grading were 
identified as rows of tiny tube-shaped cell aggregations. Trabecular structures were 





100 urn 



Figure 2 | Immunohistochemical staining for EMA in alveolar structures 
and papillary clusters of tumor cells, (a), (b): two cases with IC NST used 
in laser microdissection are shown. Alveolar structures of IC NST display 
the whole cytoplasmatic membrane and/or cytoplasm expression of EMA 
(bar 100 microns), (c): one case with micropapillary foci in IC NST is 
presented to illustrate EMA expression the stromal-basal surface of 
papillary clusters of tumor cells and an inversion of cell polarity (bar 100 
microns). NCL-EMA antibody (clone GP1.4, Novocastra) was used. 



formed by two or more rows of cells. Solid structures represented groups of 
hundreds of cells with different sizes and shapes. Discrete groups of tumor cells were 
detected as single cells or as groups of up to five cells. Alveolar structures with 
rounded shapes contained 10-30 cells. It is important to note that tumors from 
different patients may have different types of morphological structures. Detailed 
descriptions and images of different types of morphological structures were presented 
in our previous paper 7 . 

RNA isolation and cDNA synthesis. Total RNA was extracted from 69 samples of 
pre- and post-NAC tumor tissues using the RNeasy Mini Kit Plus DNase I digestion 
(Qiagen, Germany). Ribolock RNase inhibitor (Fermentas, Lithuania) was added to 
the isolated RNA. To assess RNA integrity, RIN was measured using 2200 
TapeStation Instrument and R6K ScreenTape (Agilent Technologies, Inc., Santa 
Clara, USA). RNA with RIN > 6 was reverse transcribed to cDNA using the 



Table 6 | The relationship between changes in MDR gene expression and presence/ absence of different types of morphological structures in 
breast tumors 



nl/r.2 (%/%) 







ABCB1 


ABCC1 


ABCC5 


ABCG1 


ABCG2 


Alveolar structures 


no 


7/12 (37/63) 


9/10 (47/53) 


1 0/9 (53/47) 


5/14(26/74) 


7/12 (37/63) 




yes 


27/23 (54/46) 


27/23 (54/46) 


25/23 (52/48) 


21/27 (44/56) 


24/26 (48/52) 


Trabecular structures 


no 


9/8 (53/47) 


10/7 (59/41) 


10/7 (59/41) 


7/9 (44/56) 


10/7 (59/41) 




yes 


25/27 (48/52) 


26/26 (50/50) 


25/25 (50/50) 


19/32 (37/63) 


21/31 (40/60) 


Tubular structures 


no 


15/13 (54/46) 


15/13 (54/46) 


13/13 (50/50) 


1 1/16 (41/59) 


12/16 (43/57) 




yes 


1 9/22 (46/54) 


21/20 (51/49) 


22/19 (54/46) 


15/25 (38/62) 


19/22 (46/54) 


Solid structures 


no 


20/1 T (65/35) 


17/14 (55/45) 


18/12 (60/40) 


1 1/20 (35/65) 


17/14 (55/45) 




yes 


14/24' (37/63) 


19/19 (50/50) 


1 7/20 (46/54) 


15/21 (42/58) 


14/24 (37/63) 


Discrete groups of tumor 


no 


14/8 (64/36) 


10/12 (45/55) 


13/9 (59/41) 


10/11 (48/52) 


15/7 2 (68/32) 
16/31 2 (34/66) 


cells 


yes 


20/27 (43/57) 


26/21 (55/45) 


22/23 (49/51) 


16/30 (35/65) 



Statistical analysis: Pearson's chi-squared test was used to detect significance of relationship between changes in MDR gene expression and the presence/absence of different types of morphological 
structures. 

p = 0.022 (Bonferroni-corrected p = 0.55); 
2 , p = 0.008 (Bonferroni-corrected p = 0.2). The Bonferroni-correction was calculated as the each p value multiplied by the number of tests (n = 25). 

Abbreviations: nl , number of patients with increase in gene expression after neoadjuvant chemotherapy; n2, number of patients with decrease in gene expression after neoadjuvant chemotherapy; no, the 
absence of any morphological structures; yes, the presence of any morphological structures. 
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RevertAid Kit with random hexanucleotide primers (Fermentas, Lithuania) following 
the manufacturer's instructions. 

Laser microdissection. Frozen tumor samples from two untreated patients with IC 
NST were used for PALM MicroBeam laser capture microdissection (Carl Zeiss, 
Germany). Alveolar, trabecular, solid structures, and discrete groups of tumor cells 
were isolated from 5 (J.m-think sections stained by hematoxylin (Dako, Denmark) 
and eosin {Dako, Denmark) (Fig. 1). In addition, the microenvironment of these 
structures was also isolated. Note that tumors of these patients did not contain tubular 
structures. The microdissected material was collected in RLT lysis buffer (RNeasy 
Plus Micro Kit, Qiagen, USA), and RNA was extracted according to the 
manufacturer's instructions. Ribolock RNase inhibitor (Fermentas, Lithuania) was 
added to the isolated RNA. RIN was measured using 2200 TapeStation Instrument 
and High Sensitivity R6K ScreenTape (Agilent Technologies, Inc., Santa Clara, USA). 
cDNA was synthesized, ligated, and amplified using a QuantiTect Whole 
Transcriptome Kit (Qiagen, USA) according to the manufacturer's instructions. 

Expression analysis. The expression levels of the MDR genes were measured by 
quantitative real-time PCR (qRT-PCR) based on TaqMan technology using a Rotor- 
Gene-6000 instrument (Corbett Research, Australia). qRT-PCR was performed in 
triplicate reactions in a volume of 15 jil containing 250 uM dNTPs (Sibenzyme, 
Russia), 300 nM forward and reverse primers, 200 nM probe, 2.5 mM MgCl 2 , 1 X SE 
buffer (67 mM Tris-HCl pH 8.8 at 25°C, 16.6 mM (NH 4 ) 2 S0 4 , 0.01% Tween-20), 
2.5 U Hot Start Taq polymerase (Sibenzyme, Russia), and 50 ng of template cDNA. 
Samples were heated for 10 min at 95 C followed by 40 cycles of amplification for 
10 s at 95 °C and 20 s at 60 C. The primer and probe sequences of ABCB1, ABCC1, 
ABCC5, and ABCG2 were given in our previous study 15 . The primer and probe 
sequences of ABCGl were obtained from a previous paper 44 . Two internal genes - 
GAPDH (in case of RNA from tumor bulk) and ACTB1 (RNA from the 
microdissected material) were used to normalize expression levels of the studied 
genes. The average C t (cycle threshold) was estimated for both the gene of interest, 
GAPDH, and ACTB1. Relative expression was evaluated using the Pfaffl method 45 , 
and the formula was used to determine the expression ratio between the sample and 
the calibrator 15 . The relative expression level was also normalized to a calibrator 
consisting of a pool of normal breast tissue specimens. For this purpose, specimens of 
adjacent normal breast tissue from 10 breast cancer patients (NAC-free) were used as 
a source of normal RNA. In case of the microdissected samples, we normalized the 
expression levels relative to normal breast tissue of the same patient. The results were 
presented as n-fold differences in MDR gene expression relative to GAPDH/ACTB1 
and normal breast tissue. 

Statistics. Statistical analyses were performed using STATISTICA 8.0 software 
(StatSoft, Tulsa, OK, USA). The arithmetic mean value and standard error were 
calculated for each sample group. Logistic regression was applied to identify the link 
between pre-/post-NAC expression levels of MDR genes and the presence/absence of 
different types of morphological structures in breast tumors. Pearson's chi-square test 
was used to detect the association of NAC response and change (increase or decrease) 
in expression of MDR genes with the presence/absence of different types of 
morphological structures in breast tumors. The Bonferroni correction was applied to 
address the problem of multiple comparisons and was calculated as the each p value 
multiplied by the number of comparisons 46 . The necessary to apply Bonferroni 
correction to the comparisons made was dictated by the confirmatory nature of this 
study with a clear prespecified key question consisting of several hypotheses analysed 
by multiple significance tests 47 . P-values that were corrected to values more than 1 
were truncated to 1 . Differences were significant if the corrected p value was less than 

0. 05. All p values were two-sided. 
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